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We describe a numerical technique to compute the equilibrium measure, in logarithmic 
potential theory, living on the attractor of Iterated Function Systems composed of one- 
■ dimensional afiine maps. This measure is obtained as the limit of a sequence of equilibrium 

measures on finite unions of intervals. Although these latter are known analytically, their 
computation requires the evaluation of a number of integrals and the solution of a non-linear 
set of equations. We unveil the potential numerical dangers hiding in these problems and 
j /^ | we propose detailed solutions to all of them. Convergence of the procedure is illustrated in 

^ . specific examples and is gauged by computing the electrostatic potential. 
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Suppose that a finite amount of "charge" is placed on a set of "conductors" of arbitrary shape, 
placed at equally arbitrary positions in space, and connected by thin, conducting wires. Then, this 
charge distributes itself so to minimize repulsion, i.e., the electrostatic energy of its configuration. 
It is very easy for us to intuitively perceive this phenomenon, and indeed we have been exposed 
to it since early in school, via easily realizable laboratory experiments that go back to the years 
of Alessandro Volta and to his electrophorus 1 , which he did not invent, but he perfected into 
an instrument to accumulate and transfer electric charge, years before his monumental Pila. This 
intuition is common knowledge among scientists of all disciplines: in Italy, at least until recently, 
electrostatics was taught in high school, both scientific and of classical kind and, even more 
amazingly, it was included in the syllabus of the degree in mathematics and of medical school 
alike! This paper is dedicated to a wonderful teacher who let middle-school students, twelve 
years-old, play with electrostatic devices and inspired them to pursue scientific research. 

What is far less known is the fact that the mathematical scheme behind this physical de- 
scription can be easily generalized to spaces and conductors of arbitrary dimension — even fractal 
PT?5l [Tfj] — and to arbitrary analytical forms of the "electrostatic" interaction, that need not be 
"Coulomb" , but can take on different analytical forms. These generalizations have proven to be 
extremely fruitful in a wide variety of fields, from PDE's to harmonic analysis, to name just a few. 
In this paper, we will put ourselves in the situation where the ambient space is two-dimensional 
(the complex plane) and the pseudo-Coulomb law depends on the inverse of the distance be- 
tween two charges: physicists will easily think of the interaction of two infinite, parallel wires, 
mathematicians will immediately recognize the logarithmic potential that takes us into the realm 
of the well established logarithmic potential theory UOj . 
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Indeed, if we let E to be a compact domain in the complex plane C, and a a positive 
probability measure supported on E, the potential V(a; z) generated by a at the point z in C, 
is given by the formula 

V(a;z) :=- [ \og\z-s\ da(s). (1) 

JE 

The electrostatic energy £{o~) of the distribution a is then given by the integral 

8(a) := I V(a;u)da(u) = - f [ log \u - s\ da(s)da(u). (2) 
Je Je Je 

The equilibrium measure <je associated with the compact domain E is the unique measure that 
minimizes the energy £ (er), when this latter is not identically infinite. Potential theory gives us 
a formidable set of technical instruments [7j [TU] to deal with this problem in relation to general 
compact sets E. 

In this paper, we are interested in a specific family of compact sets E, that are Cantor sets on 
the real line associated with the construction of Iterated Function Systems (IFS) — to be defined 
momentarily — and we study the efficient numerical computation of their equilibrium measure <je 
and of its potential V(a; z). Our results will consist of a reliable algorithms to achieve these goals. 
This algorithm is rather straightforward, in the sense that it combines well known theoretical 
facts (the analytical form of the equilibrium measure when E is a collection of intervals) and well 
known numerical techniques (Gaussian summation, root finding) with a standard idea in IFS 
construction: the hierarchical approach. Despite or rather because of this simplicity, numerical 
experiments indicate that the algorithm is stable and it can be used as a microscope to probe 
deeply into the structure of the equilibrium measure on Cantor sets. 

This paper is organized as follows: in the next section we define the attractors of Iterated 
Function Systems, via a hierarchical construction. Each step in this construction yields a finite 
family of intervals, whose equilibrium measure is analytically known. We sketch this solution in 
Sect. [3] The following three sections contain the numerical techniques to find this solution and 
to overcome a few potentially destructive numerical difficulties. In Sect. 0]we apply Gaussian 
integration to compute a class of integrals that reduce the problem to a set of non-linear equa- 
tions, and in Sect. [5]we describe a convenient approach for the solution of these equations. Next, 
in Sect. [6j we show how to tailor the root finding routine to the hierarchical structure of IFS 
attractors. As a result, we obtain the equilibrium measure on these latter. Pictures and graphs of 
numerical experiments illustrate the theory. In section [7] we compute the potential generated in 
the complex plane by the equilibrium measure on the IFS attractor and we compute its capacity. 
The application of this approach to deep questions on equilibrium measures on Cantor sets is 
briefly discussed in the conclusions. 



2 Attractors of Iterated Function Systems 

Let us therefore construct the set E of which we want to place an electrostatic charge. To do 
this, we need to recall, as briefly as possible and in the simplest setting, the construction of 
Iterated Function Systems (IFS) [26j [T9l QH H2 [13] . These are collections of maps 0, : R -> R, 
i = 1,...,M, for which there exists a set A, called the attractor of the IFS, that solves the 
equation 

A = |J <f>i{A):=*(A). (3) 

j=t,...,M 

When the maps (f>i are contractive, the attractor A is unique. It can also be seen as the fixed 
point of the operator defined in eq. ([3]), on the set of compact subsets of R. Since this space is 
complete in the Hausdorff metric, and since $ turns out to be contractive in this metric, the set 
A can be also found as the limit of the sequence & n (E°), where E° is any non-empty compact 
set: 

A = lim $ n (E°). (4) 

n— too 
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Attractors of Iterated Function Systems feature a rich variety of topological structures, so 
that their full characterization is far from being fully understood, especially in the case of IFS 
with uncountably many maps [21] [22j [17j [18] . We nonetheless restrict ourselves in this paper to 
a particular family, that of IFS with a finite number of affine maps of the form: 

<^( s ) = 6 j( s - 7j) + 7i> 3 = h ■ ■ ■ ,M, (5) 

where Si are real numbers between zero and one, called contraction ratios, and ji are real con- 
stants, that geometrically correspond to the fixed points of the maps. Under these conditions, 
the attractor A is a finite or infinite collection of intervals, or a Cantor set. We will consider 
this last, interesting case. It can be easily recognized that the famous "middle-third Cantor set" 
follows in this class: it can be obtained by taking just two maps with 6 — 1/3 and 71 = — 1, 
72 = 1. We shall study this case later in this paper. 

As remarked above, any non-empty compact set E° can be used in eq. ((3]). We let E° be the 
convex hull of the attractor A, that can be easily be identified as the interval 

£° = [7i,7m], (6) 

where we have ordered the IFS maps according to increasing values of their fixed points: jj < 
7j+i, for any j = 1, . . . , M — 1. This is the core of the common hierarchical construction of the 
set A, as the limit of the sequence of compact sets E n : 

M n 

E n = $ n (E°) = \J[a i ,p i \. (7) 

i=l 

In the case of fully disconnected IFS (i.e. those for which the intervals 4>j(E°) are pairwise 
disjoint), E n is the union of N := M n disjoint intervals, that we will denote as _E™ := [c*j,/3j]. 
In other words, the set E n can be seen as a vector of intervals, stored via the two vectors of 
their extreme points. Of course, these latter should also carry a superscript n, referring to 
the generation. Not to overburden the notation, we will leave this superscript implicit, when 
confusion is not possible. We restrict ourselves to the disconnected IFS case, for the remainder 
of the paper. 

For similarity with the spectral analysis of periodic solids, we will call the intervals Ef at 
r.h.s. of eq. bands at generation (or level) n. It appears from eq. (0 that all bands Ef 
at level n can be obtained by applying the transformations <f)j, j — 1, . . . , M, to the bands at 
generation n — 1. As a consequence, the length of these intervals, as well as their sum, tend to 
zero geometrically and the attractor is a Cantor set of null Lebesgue measure. 

While none of the bands -E™ -1 are also bands at level n, the contrary happens for the so-called 
gaps between the bands. For these, we will use the notation: 

M n -i 

G n :=\jG?= |J (A,a i+ i). (8) 

i i—1 

At level zero we have no bands, while at level one these are in the number of M — 1. It is 
then noticeable that all gaps at level n — 1 remain gaps at level n, and at the same time a set 
of new gaps is generated from the former: let H n be the set of "new" gaps created at level n, so 
that 

G n = G n - 1 \jH n . (9) 
The latter set, H n , is iteratively constructed as 

M 

/f«=U^(if"- 1 ) ) (10) 

j=l 

starting from H 1 = {Px,ot\) [_}■■■ U(/^m-i> a \i)- This seemingly complicated algorithm is noth- 
ing more than the straightforward generalization of the construction of the ternary Cantor set 
by deleting the "middle third" . This property of the set of gaps is important in our computation 
of the equilibrium measure on the Cantor set A. 



3 



3 Equilibrium measure on a set of intervals 



Let us now embed E n in the complex plane, and study its equilibrium measure. The solution 
of the equilibrium problem for a finite union of N intervals [aj, i = 1, . . . , N, is well known 
HO [251 HI- Define the polynomial Y(z), 



N 



Y(z)=Y[(z-a l )(z~(3 i ), (11) 



and its square root, y/Y(z), as the one which takes real values for z real and large (larger than jm 
indeed). Also, let the real number Q belong to the open interval (ft, «i+i), for i = 1, . . . , N — 1 
(i.e. to the gap Gf). Define Z(z) as the monic polynomial of degree N — 1 with roots at all Q's: 



N-l 



Z(z)=H(z-C i ). (12) 

i=l 

With these premises, there exists a unique set of values i = 1, . . . , N — 1} that solve the set 
of coupled, non-linear equations 

Oi+l Z(s) 

. w ds = 0, i = l,...,iV-l. (13) 

Our first task in the following will be to evaluate numerically these integrals for any given set of 
values {Ci,i = 1, • ■ • , N — 1}. The second task will be to find the solution of this set of equations. 
This is of paramount importance, for it permits us to find the equilibrium measure of the set E n : 
for simplicity of notation, we denote it by a n := oe^ '■ 

d ^^) = ^X>NA](«) da. (14) 



i=l 



It is apparent from the previous equation that the measure a n is absolutely continuous with 
respect to the Lebesque measure on E n . Of particular relevance are also the integrals of <r" over 
the intervals composing E n , that we call the harmonic frequencies: 

:= ~ I"' -PL*- (16) 

In fact, when these frequencies are rational numbers, of the kind Pi/N, with pi integer, there 
exists a strict-T polynomial on E n [3], that is, a polynomial with oscillation properties mimicking, 
and extending, those of the classical Chebyshev polynomials on a single interval. Also, we will 
consider the integrated measures ft" n : 



nS,:=5>?, (16) 



i=l 



that are the integral of a n on the interval [71, f3 m }. 

Finally, the Green function with pole at infinity can be computed as the complex integral 

"M-jCt^T* (17) 

in which no absolute value appears. The Green function can be used to compute the logarithmic 
capacity of the set E n , C(E n ), via a real integral 

los(C<JS " ))= £ [ ^j-^m)i ds ' (18) 

Although our technique could permit to evaluate the above integral, we will compute the capacity 
C(E n ) following a different approach. 
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4 Integrating the ratio Z/y/Y 



The first step in the determination of the equilibrium measure is the numerical determination of 
integrals of the kind (TT3"|) . or of the kind (|T5|) . that are quite similar. A straightforward technique 
to achieve this goal is described in this section. The first integrals are: 

£,:=!/" -JM= da, (19) 

where Gi is a gap at generation n. The case when Gi is replaced by a band Ei can be handled 
by minor variations of the technique that we are going to describe. For simplicity of notation, 
since n is fixed, we will omit its mention in the following derivation. A full set of integrals /Cj, 
for variable i, need to be evaluated. Let us now keep the index i fixed and show how to compute 
any single one of them. 

Firstly, consider the change of variables induced by the affine transformation tpi :— x — > 
AiX + Bi, that maps the interval Gi into E° = [71, 7m]- Without any loss of generality we can 
assume that E° = [—1,1]. The coefficients Ai and Bi can be computed by careful bookkeeping 
from the map parameters Sj and /3j. Even if we feel no need to report the explicit formulae 
here, they are essential in the numerical implementation that we will describe in the following. 
In so doing, the bands E m and the gaps G m , for m different than i, are mapped, via ipi, to 
intervals outside [—1, 1]. Equivalently, the variables £ m are mapped via the same map to new 
values tpi(Cm)- New functions Zi and Yj, are then written as in eqs. (jlip and (|12[) . now in terms 
of the transformed V'i(Cfc), ipi{phn) an d ipi(flm)- In conclusion, the integral JCi can be written as 

Secondly, remark that the product 1 — z 2 appears in the rescaled function Yi(z). Therefore, 
we can part the full product Yi(z) in two factors: 

Y l (z) = (l-z 2 )Y l (z), (21) 

where Yi(z) is implicitly defined. With these notations, it appears clearly that the integration 
of f'J s ^ with respect to the Lebesgue measure on f— 1, 11 can be seen as an integration with 

respect to the Chebyshev measure ds/ivyl — s 2 . Then, the integral K% can be approximated by 
a Gaussian summation, 

f 1 Zj(s) ds Zj{x k ) ( . 

where x k , Wk are Gaussian points and weights for the Chebyshev measure, respectively. As a 
consequence, each integral /Q can be readily evaluated by Gaussian summation. 

Clearly, the above assumes that the values Q be known for alH = 1, . . . , N — 1. To compute 
these latter, the N equations K-i = 0, i = 1, . . . , N — 1, must be solved. In the next section 
we employ an algorithm for the solution of this set of non-linear equation that requires the 
computation of the derivatives dJCi/dC, m . Also these quantities can be computed via Gaussian 
summation. Observe in fact that the partial derivative of eq. (|22p requires the computation of 

1 ' 

Z l (z) = -A l Y[(z-^(Q)), (23) 



that therefore leads to 



did . A T\i^ m (xk - Mb)) .... 

— ~-Ai^w k . (24) 

dQm *=i J\Y(x k )\ 



5 



5 Solving the non-linear equations /Q = 



We have seen that the determination of the equilibrium measure on the set E n requires the 
solution of a system of N — 1 non-linear functions in N — 1 variables, /C^ (Ci , ■ - ■ , CiV-i) = 0. 
We can compute this solution by a careful usage of the routine HYBRJ in Minpack 2 , which 
employs a modification of the Powell hybrid method. This technique is a variation of Newton's 
method, that iteratively converges to the solution vector. As it turns out, HYBRJ cannot 
be applied blindly without taking into account the particular nature of the set E n and of the 
discretized integrals (|2"2"1) . In fact, in the construction of the attractor of a fully disconnected 
IFS, all bands E™ have a length that goes to zero geometrically fast, while gaps differ in length 
by orders of magnitude. This requires a series of procedures, in the evaluation of the integrals 
K,i and that we explain in this section. 

The first consists in the approach described at the end of the preceding section: it is convenient 
to scale the each interval G™ to [—1,1], when performing the relative integration, to employ the 
Gaussian form in eqs. (|2"2")l and ([2"1|) . 

When this is done, a further difficulty is to be overcome: the function Zi(xk) / \J \Yi(xk)\ is 
made of factors differing by various orders of magnitude, that must be conveniently rearranged 
when evaluating it numerically. Therefore, in the expansion 

\J\Yi{x k )\ \J\ X k ~ ^{ai)\\Xk ~ MPrn+l)\ Um^i,i+1 \ X k ~ M^m)\\x k ~ MPrn)\ 

we find it convenient to group together the ratios (xk — V'i(Cm))/ vf^fe ~ ipi(&m)\\xk — ^i{Pm)\, 
for m < i, and (x k - ipi(( m ))/ V\ x k ~ ^i{a m +i)\\xk ~ ipi(fi m+ i)\, for m > i, in which numerator 
and denominator have comparable size. The more the index m is different from i, the larger are 
the differences between xu and the points ipi(a m ) and ipi(/3 m ), and the more these ratios tend to 
one. Physically, this means that remote gaps GJ^ have a small influence on the integral Ki. our 
ordering correctly reproduces this fact. As it appears from eq. (125p . grouping terms as above 
leaves out a few factors, those close to the interval Vi(G") = [—1,1], that mostly contribute to 
the value of the full product, and are therefore taken care separately. A similar procedure can 
also be applied to the formula for derivatives, eq. (pM)) . 

But the more important precaution to be taken follows from considering the variables involved 
in the Newton- like iteration of HYBRJ. In the original equations K4 = the unknowns are the 
positions of the N—l zeros of Z, i.e. f m , m = 1, . . . , JV — 1. Each of these belongs to the gap G™ . 
It is then apparent, from the description in Sect. [21 that the range of variation of these values 
can differ by orders of magnitude, and can quickly become smaller than numerical precision. It 
is then mandatory to resort to new variables: we consider in place of Cm the normalized variables 
A m := VVra(Cm): where, as above, ip m is the affine transformation that maps G^ into [—1,1]. 
These new variables are all bound to the interval [—1,1]. Some care must now be taken in 
the evaluation of the compounded derivatives, but fortunately, the new variables are the most 
natural, in the construction of the IFS, as it can be easily recognized from eqs. ©,([71). 

The above provides us with a working code for the computation of the non-linear functions 
(|22p and of the partial derivatives ([2"B"|). as functions of the variables A TO , m — 1, . . . ,N — 1, that 
can be linked as a subroutine to the fortran (yes, fortran!) code HYBRJ, run in double precision 
on a 32-bit processor. In Figure (fT]) we show the absolute values of the integrals Id before and 
after the call to HYBRJ in a typical case, that of the ternary Cantor set at generation number 
n = 7, to be defined more precisely later on in Sect. [JJ The algorithm has converged to the 
correct solution within numerical precision. 

To fully appreciate the precision of the solution vector, we can examine the value of the 
partial derivatives S^- at the solution vector. These values are reported in Figure ([2]) versus the 
difference i — m: as remarked above, the larger this difference (in absolute value) the lesser the 
influence of X m on the integral A^: this is clearly evident in the figure. Also, we can observe that 
the system of equations is almost "diagonal" , in the sense that diagonal components with i = m 
(the tip of the "Christmas tree") are orders of magnitude larger than the non-diagonal branches, 
i ^ m. 
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Figure 1: Absolute value ti of the non-linear functions K-i at the initial vector (Ax, . . . , Ajv-i) (see 
next section for its definition) (green crosses), and after the action of the root finding routine 
(red crosses). The case under study is the ternary Cantor set at generation number n = 7. The 
number of Gaussian points employed is 2048. 
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Figure 3: Roots A m of the equations (fT3|) versus generation number n, for m = 1, . . . , 2™ _1 — 1, 
for the IFS described in the text. Lines connect roots in the same gap at different n. The number 
of Gaussian points employed is 160. 

6 Hierarchical structure of the equilibrium measure 

We have described in the previous section an algorithm for the determination of the roots of the 
equilibrium equations (I13[) . The algorithm applies to any finite family of intervals. In particular, 
we want to apply it to a system of intervals generated by a fully disconnected IFS. In the previous 
section we have already presented results in Figures [1] and [2] referring to the case of the middle- 
third Cantor set. A further remark on the way they have been obtained is important. While the 
case of the ternary Cantor set is particularly favorable because of its symmetry, we now choose 
to describe our procedure in the case of an asymmetric IFS given by two-maps with 6% = 4/5 
and 62 = 1/10. We will return to the middle-third Cantor set in the next section. 

Observe that the cardinality of E n is M n , where M is the number of IFS maps, and therefore 
the number of gaps is M n — 1. Because of eq. ©, a subset of cardinality M "~ 1 — 1 of these latter 
exist already as gaps at generation n — 1. The old gaps are approximately (M — 1) times less 
numerous than the new ones, nonetheless, their role is crucial. In fact, the Newton-like search of 
the algorithm HYBRJ in Minpack need to be initialized. It is particularly efficient to initialize 
the values of the roots A^ /m to the "old" values at generation A^~ , and the new ones to zero 
(which means that ( m lies in the middle of G m .) 

Everything described is clearly evident in Fig. [3l that plots the values of the solution variables 
A m versus n. We can observe the new gaps that are created at each new generation: data are 
drawn as crosses, and lines join the A values of roots of Z lying in the same gap for successive 
values of n. As a gap gets "old" , the root inside it tends to a limit value. From Fig. [3] we can 
also observe that all roots have absolute value less than 1/10, while in principle they are bound 
by one. 

After all this preparatory work we are now ready to examine the equilibrium measure on the 
set E n . It is straightforward to obtain the harmonic frequencies and the integrated measures 
O^i, discussed in the second section, via the Gaussian integration technique described in Sect. 
[4] The case of fig. [3] is examined again in figure [4] by plotting the integrated measures OJ^ versus 
n: recall that fij^ is the measure of the set (71, Cm) under the equilibrium measure on E n . As 
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Figure 4: Integrated measures f2J^ versus versus generation number n, for the IFS described in 
the text. Lines connect points as explained in the text, and as done in the previous figure [3J 
They are only apparently horizontal (see the following Fig. |6]) . The number of Gaussian points 
employed is 160. 

such, these values are the heights of the plateaus in Fig. [5j that draws the integral 



versus x, for the same IFS described above. 

As n grows, E n tends in Hausdorff distance to the attractor of the IFS, that is a Cantor set 
of zero Lebesgue measure. We want now to show that our technique can be used to compute 
the equilibrium measure of the IFS attractor itself: it appears from Fig. [5] that 0, n (x) tends to a 
limit function Q(x) for Lebesgue a.e. x. This function is precisely the integral of the equilibrium 
measure on the IFS attractor. 

In fact, we have already remarked that gaps, once created, stay forever in the complementary 
of the Cantor set. Let us therefore gauge the convergence, as n grows, of the equilibrium measure 
of the infinite interval to the left of each gap. This is nothing else than fij),( n ), when the index of 
the gap, m, depends on the generation n according to the rule m(n + 1) = Mm[n). In Figure |6] 
we choose the first gap, G\, that is labeled at successive generations as GV^\ we extract a single 
line from fig. [3] and we replot it alone so to render evident its variation. It appears that we have 
exponential convergence of the equilibrium measures, to the limit equilibrium measure of the 
Cantor set. We have found the same behavior for all gaps, so that we can numerically observe 
convergence to the limit function Q,(x) . But measures are born for integration, and integrating 
is what we are now up to. 



The convergence that we have observed at the end of the previous section is just an instance 
of the convergence of a n to the equilibrium measure on the IFS attractor. A second example 
is obtained considering the electrostatic potential V(a n ;z) defined in eq. (fT]). Our technique 
permits to compute it in the full complex plane. Care must be exerted when z belongs to E n , 




7i 



(26) 



7 The Capacity of E n and of the IFS attractor 
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Figure 5: Harmonic measures tt n (x) versus x at various generation numbers n for the IFS 
described in the text. The number of Gaussian points employed is 160. 
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Figure 6: Integrated measure ft^n versus generation number n for the IFS described in the text 
(red crosses). Also plotted is the fit of the discrete values by the curve f(n) = a + be~ cn (green 
line). The number of Gaussian points employed is 160. 
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Figure 7: Electrostatic potential V(a n ; x) versus x, at generation level n = 1, . . . , 5, for the IFS 
generating the ternary Cantor set. The number of Gaussian points employed is 2048. 

for it may coincide with one of the Gaussian points, in which case one simply has to change their 
number in order to remove the coincidence. In this section we consider again the IFS generating 
the middle-third Cantor set, given by two maps with S = 1/3 and 71 = — 1, 72 = 1. 

In Fig. [7]we plot the potential V(a n ; x) when x belongs to the interval [0, 1], for n = 1, . . . , 5. 
Notice that for symmetry we have constructed our ternary Cantor set so that its convex hull 
is [—1,1]. In the range of the figure, we compute the potential both on the set E n and on its 
complement. In the gaps, as the generation level n grows, we observe exponential convergence to 
a limit value. At fixed n, the potential V(a n ;x) must take a constant value almost everywhere 
on E n . This is also observed in Fig. [7] This constant value is linked to C(E n ), the capacity of 
E n , by 

V(a n ; x G E n ) = - log(C(£™)), a n a.e. (27) 

Our technique therefore also permits to compute these capacities, by computing the potential. 

In Table 1 we plot V(a n ; x) versus n for a value of x in E n , n = 7, close to 71 (second column), 
together with the average potential, J2i V(cr n ; xi)/L, taken over L points in E 7 (third column). 
Observe that since E n C E m for any n > m, this choice guarantees that the same set of sample 
points can be used for all values of n in the Table. Difference between the two columns is due, 
and can serve to gauge, both the coarseness of Gaussian integration and the numerical errors 
involved in our technique. We observe a discrepancy that steadily diminishes as n increases, from 
2 x 10~ 4 for n = 1 to 3 x 10~ 6 at n = 7, in the same conditions. 

In Fig. [7] we observe convergence of the potential V(a n \ x) for increasing n at fixed x, also for 
x G E n , although at a visibly slower rate than for x in the complement of E n , This is particularly 
interesting if we now try to extrapolate the values in Table 1 to evaluate the capacity of the 
attractor of the IFS: Figure [8] displays the same values of the second column of Table 1 together 
with a fit by an exponential function. Data from the third column would be indistinguishable 
from these, on the scale of the figure. Although much better extrapolation algorithms can be 
applied [23], this simple technique already provides an asymptotic value of the potential, which 
implies, via eq. (|27p . an approximate value of the capacity C(A) = 0.44189238 that agrees within 
6 x 10~ 6 with the value 0.441898204379014 provided by [6] (obviously doubled, because of the 
different scaling employed here). Performing the same extrapolation technique on the average 
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n 


K(o- n ;a:) 


(V(a n ;x)} 


1 


U. 751045 


U. 752051 


2 


0.799586 


0.783380 


3 


0.807899 


0.807941 


4 


0.807899 


0.807941 


5 


0.812184 


0.812210 


6 


0.814388 


0.814392 


7 


0.815506 


0.815509 



Table 1: Potential values V(a n ;x) at x = -0.999996236647154 and average values (V(a n ;x)) 
over L = 4096 points for the middle-third Cantor set IFS at different generation levels n. The 
number of Gaussian points employed is 2048. 




0.75 - 
0.74 - 
0.73 - 



0.72 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 7 

n 

Figure 8: Data from the second column in Table I: potential values V(a n ;x) at x = 
—0.999996236647154 versus n (red crosses). Also plotted are the fit of the discrete values 
from n = 4 to n = 7 by the curve f(n) = a + be~ cn (green line) and the infinite limit value 
a = 0.81668890 (blue line). The other fitting parameters are b = -0.1278376 and c = 0.66927525. 
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values in the third column of Table 1 yields an even better estimate, C(A) — 0.44189726, this 
time within 1 x 10~ 6 of the exact value. We can therefore conclude that we can reliably compute 
integrals with respect to the equilibrium measure on a Cantor set. 

8 Conclusions 

We have presented in this paper a technique for computing the equilibrium measure on attractors 
of iterated function systems. This technique consists of the careful concatenation of different 
ideas in numerical analysis, potential theory, and IFS construction. We have tested our approach 
by computing the capacity of the middle-third Cantor set, for which a reliable independent 
computation exists. 

Yet, our interest in this problem does not end here: we intend to use the techniques of 
these paper as a mathematical electrophorus to compute the Jacobi matrix of the equilibrium 
measure. In fact, since this matrix embodies the orthogonal polynomials of the measure, its 
asymptotic properties are of paramount importance in a number of problems of importance in 
harmonic analysis, quantum mechanics and fractal geometry [HI HH EH EH EH El EH]. We 
have reasonable hope that this investigation will help us to clarify a long-standing conjecture on 
the almost-periodicity of Jacobi matrices associated with IFS measures [50] . 
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